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Abstract. Motivated by various spin-1/2 compounds like CS2CUCI4 or k-(BEDT- 
TTF)2Cri2(CN)3, we derive a Raman-scattering operator a la Shastry and Shraiman 
for various geometries. For T=0, the exact spectra is computed by Lanczos algorithm 
for finite-size clusters. We perform a systematic investigation as a function of J2/J1, 
the exchange constant ratio: ranging from J2 = 0, the well known square-lattice case, 
to J2/J1 = 1 the isotropic triangular lattice. We discuss the polarization dependence 
of the spectra and show how it can be used to detect precursors of the instabilities of 
the ground state against quantum fluctuations. 



1. Introduction 

Highly frustrated magnetic systems are highly susceptible to quantum spin fluctuations 
and instabilities towards competing ground states. The triangular Heisenberg 
antiferromagnet with spin S = 1/2 constitutes a paradigm of that class of systems. 
In that context, it is interesting that experimental investigations of the CS2CUCI4 pQ 
and /t-(BEDT-TTF)2Cu2(CN)3 [2j materials, which can be both described as a first 
approximation by a triangular lattice with spatially anisotropic exchange couplings, 
indicate exotic behaviors. Exotic behaviors include either the realization of a spin-liquid 
ground state or a magnetically ordered phase with a magnetic excitation dispersion 
strongly renormalized compared to the classical spin-wave. Indeed, recent numerical 
studies based on series expansion [2] have found that frustration manifests itself rather 
directly in the spin excitations of archetype models of two-dimensional frustrated spin 
systems. A softening of the magnon excitations in a broad region of the reciprocal q 
space is observed as a quantity, /, which parametrize the level of frustration, is increased. 
Since, within a semi-classical picture, the effect of frustration is to bring about the 
effect of competing ground states, it is expected that the spin excitations out of a semi- 
classical long-range ordered ground state are direct tell-tale indicators of the presence 
of frustrating interactions. In this paper we explore the possibility that frustration can 
also be investigated via the two-magnon density of states at zero wavevector measured 
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by polarized inelastic magnetic Raman scattering. 

The rest of this paper is organized as follows: in the first Section we present the 
derivation of the scattering operator, then we discuss its form for various polarizations. 
In the last part we present the associated Raman spectra obtained by exact- 
diagonalization of finite-size clusters. 

2. Scattering operator 

We first consider the Hubbard model on the anisotropic triangular lattice: 

H = Hk + Hu = ^2 *w ( C L C J> + h - c -) + U n ^ Ui ^ ( 1 ) 

where is the hopping from a site i to a neighboring site j, a is the spin degree of 
freedom, and c and c' are the usual annihilation and creation operators. For the case 
considered here, there are two distinct hopping integrals, t\ and t2, along the corre- 
sponding directions shown in Fig. ^ 




Figure 1. Geometry of the lattice and polarization vectors. Left: The dashed lines 
represent the coupling J2, the thick lines are for the exchange coupling J\. Right: 
Definition of the polarization vectors by two angles 4> and 6. 



At half-filling and in the large-t/ limit, the system is described (to second order in 
t/U) by an effective spatially anisotropic Heisenberg Hamiltonian: 

7i c s = JijSi ■ Sj (2) 

where the different exchange coupling along each direction is due to the difference in 
the hopping amplitude, i.e J\ = At\/U , J2 = &t\/U. 

Raman scattering consists of an incoming photon (of energy uji) scattered into an 
outgoing photon (of energy ujf), involving different manifold of electronic states having 
zero or one double-occupancy. These transitions depend on the polarizations (referred 
as Ej n and E out ) of the incoming and outgoing photons. Following the early work of 
Fleury and Loudon jl] and Shastry and Shraiman [S], we derive an effective scattering 
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spin Hamiltonian describing this problem. The a-component of the electronic hopping 
current operator is: 

j« (q) = i (f^) eiQ ' (r+f ] K(r + ")ca(T) - 4 (r)c CT (r + „)] (3) 

r,i/=±(ei,e2,e3) 

where q = k/ — kj is the momentum transfer, and the band-energy. In what follows 
we take the most general set of polarization vectors, as shown in the right panel of 
Fig. HJ 

cos0 x + sin0 y 
cos 9 x + sin 9 y 

The Raman scattering process involves an energy transfer Q — Uf — Ui, the momentum 
transfer being set to zero. Therefore, we restrict ourselves to the assumption k/ « kj « 
0, and we have: 

j • E in = % cos (f> [-ti (4c r+ei - h.c) -| (4c r+e2 - h.c) 

+ y (4 C r+e 3 - ^-c)] 
Sin0 [t 2 (4 C r+e 2 - ^- c ) + h (4 C r+e 3 ~ ^-c)] 

(5) 

j ■ E out =icos9 [-tt (4c r + ei - h.c) -f (4c r +e 2 - fr-c) 

+ f (4 C r+e 3 - h.c)] 

-i^Y sin 61 [t 2 (4 c r+e 2 - h.c) + t x (c\c r+e:i - h.c)] 

The Raman scattering operator is given by the second order formula, \i) and |/) being 
respectively the initial and final state, being an intermediate state: 

Following the same algebra steps as in [Hj, and restricting \i) and |/) to the manifold 
of singly occupied states, and intermediate states \fj) to the manifold of one double 
occupancy, we use the identity, | — S» - Sj = ^2 aa i \ c i,a c j,a c j !(T > c i,<j'j an d obtain the 
scattering operator in terms of spin operators. We note that within our approach the 
total scattering operator only contains terms of the form: 

a cc |s, • S 4+e „ (7) 

From Eq. JBJ) we find that the scattering Hamiltonian prefactors which will be in 
front of Eq. (jZJ) depend on the polarization vectors orientation. These prefactors 
are proportional to the exchange coupling J, therefore, in its general form, the total 
scattering operator will depend on both J\ and J 2 . 

3. Polarizations 



The two angles (p and 6 with respect to the x-axis (see Fig. |TJ) define the polarizations 
involved in the scattering process. Therefore the scattering operator depends on a 
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projector V V {Q,<$) that defines the polarization set-up. The Raman operator takes the 
form: 

H L f(0, 0) OC J^i {^l S * ' S *+ei COS 9 COS 

+ J 2 Si ■ S i+62 [cos(0 + 0) + >/3 sin(0 + 0) + 4 sin sin 0] (8) 

+ JiSi • S i+e3 [cos(0 + 0) - \/3sin(0 + 0) + 4 sin sin 0] } 
Which can be written in the compact form: 

H LF (6, 0) ex J » V »( e > ^ ■ S *+e„ (9) 

In order to compare with the square-lattice case, we now focus on the following 
polarization geometries: 



7~Llf(-^, — f ) oc J2i Ji l^i ■ Sj +ei + Sj ■ S 

7~(-Lf{^qi f ) OC Jl [Sj ■ Sj +ei — Sj • Sj +e3 



, .. ^ , e c o o l+G3 i (10) 



When the diagonal bond J2 is 0, the first line of Eq. (jlUj) is the A lg Raman operator, 
while the second line gives the B l9 Raman operator, both on a square lattice. Most 
importantly, we note that the scattering operators depend on the ratio J\j Ji if one 
takes: 

7~Clf(^,—^-) oc ^ [JiSj • S i+ei + J 2 Sj • Sj +e2 ] 
7~(-lf{\-i ^p) oc Y2i [<ASj • Sj +ei — ^Sj • Sj +e2 ] 
We would like to attract the attention of the reader to the fact that the standard parallel 
polarization does not give a straightforward A l5 -like scattering operator as it is the case 
for the square lattice, but a linear combination of Eqs. El and ^2 



4. Results 



As a first step in exploring the Raman scattering in the Heisenberg S = 1/2 anisotropic 
triangular lattice antiferromagnet, we perform exact- diagonaliz at ions for the Raman 
operator Eq. While Raman scattering studies on triangular lattice compounds like 
cobaltites have been carried out [H] , these have focused on phonons, and there are yet no 
systematic study of polarization dependent electronic Raman spectra. Hence, our aim 
at this stage is not to obtain a quantitative description of the scattering spectra for a 
specific material, but rather identify the key qualitative features emerging in the Raman 
spectrum of the anistropic triangular lattice, possibly to motivate Raman studies and 
to connect with recent neutron data PP, NMR |2j and angle-resolved photoemission [7j. 
Motivated by the successes of exact-diagonalization to identify the essential qualitative 
features of the Raman spectrum on the square lattice [HEIj and in order to compare with 
these well-known results, we explore in this study the behaviour of a 16-site cluster as 
the frustrating J2 coupling is increased from zero. The results are summarized in Fig. 121 
The general trend observed in these plots is that a softening progressively develops as the 
system becomes more frustrated. The bottom left panel of Fig. El shows, for J2 = 0, the 
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Figure 2. Exact diagonalizations on a 16-site cluster. Raman spectra for various 
values of (9, <j>) and J2/J1. The bottom (top) left panel would correspond to A\ g {B\ g ) 
polarization on the square lattice in the case Ji = 0). 



well known zero Ai g Raman scattering on a square lattice, since the scattering operator 
commutes with the Hamiltonian [9 . As J 2 increases non-commuting contribution to 
the scattering intensity increases and the spectrum develops more structure. For the 
bottom right panel, a weak intensity if observe because the scattering operator does not 
commute with the Hamiltonian even for J 2 = 0. 

From the Hamiltonian in Eq. |21 it is obvious that the ground state will depend on 
the ratio J2/J1, ranging from a Neel-like order on the square lattice when J 2 = 0, to 
a three sublattice long-range order for J 2 = J\ [TQJ [TTJ [T2]. These different type of 
order are characterized by different magnon dispersions, therefore the Raman spectra is 
expected to depend on J2/J1. In a recent paper, Zheng et al. j3] showed, using series 
expansions, that the magnon dispersion for the anisotropic triangular lattice Heisenberg 
model exhibits a roton-like minimum. This local minimum sits at the (it, 0) point for a 
square lattice with one diagonal bond and is getting softer as the ^-diagonal frustrating 
coupling increases. In the case of the standard square lattice, for the B l9 channel with 
crossed polarizations, the electromagnetic field couples to excitations along the (%, 0) 
direction. We note from [3| that this softening is a multi-magnon process that cannot 
be captured to lowest order in a 1/S expansion. However, since we are using exact- 
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diagonalization which, by its nature, treats processes for the length scale considered 
exactly, we are able to detect indications of this softening. For specificity, consider the 
Bi 9 channel in the top left-hand panel of Fig. 2. A signature of the softening of the 
magnon excitation is clearly manifest as the frustration ratio / = J2/J1 is increased 
from zero, with a shift of the spectral weight, and the main peak moving from u ~ 3 J\ 
when J 2 = to to ~ O.6J1 when J 2 = O.8J1. 

5. Conclusions and perspectives 

The results presented here show how frustration can dramatically alter the Raman 
spectrum of an otherwise non-frustrated system. In particular, we observe a frustration- 
driven spectral downshift, the landmark feature of frustrated systems. This shift is 
another indicator of the frustration-induced softening of the magnon dispersion recently 
predicted by others. From a strictly theoretical point of view, the analysis above sets the 
stage for a more complete approach. Some exact-diagonalizations on larger cluster and 
a spin-wave approach including magnon-magnon interactions, along the lines pursued 
for the square lattice JHJE], are currently being carried out. 
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